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Abstract: 



p Using short distance QCD methods based on the operator product expansion, we 

p • calculate the J/ ip photoproduction cross section in terms of the gluon distribution function 

D . of the nucleon. Comparing the result with data, we show that experimental behaviour 

^ ' of the cross section correctly reflects the x-dependence of the gluon distribution obtained 



from deep inelastic scattering. 



The charmonium ground state J/ip is much smaller than the conventional hadrons 
constructed from light u and d quarks, and much more tightly bound: we recall that 
rj/^ ~ 0.2 fm <C ^^qcd and 2Md - Mj/^ ~ 0.64 GeV > Aqcd- As a result, when a 
J/ip interacts with a 'light' hadron, it is expected to probe the local partonic structure 
of the latter, not its size or mass. The aim of this paper is to show that measurements 
of J/'?/'-photoproduction on nucleons confirm this expectation and give a rather precise 
reflection of the gluon distribution function of a nucleon. 

Using vector meson dominance (VMD), we shall first relate forward J/?/'-photopro- 
duction to J/'?/'-nucleon scattering. Next, the total J / ip-'nudeon cross section is expressed 
in terms of nucleonic gluon distribution functions, making use of the short distance QCD 
methods based on sum rules derived from the operator product expansion (OPE). The last 
missing element, the ratio of real and imaginary J/ip—N scattering amplitudes, is obtained 
by a dispersion relation. We can then study directly the interrelation of the energy 
dependence of J/V'-photoproduction and the x dependence of the gluon distribution in 
the nucleon. 
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The use of the J/tp as a probe of the partonic status of a given medium is of particular 
interest for the study of colour deconfinement in nucleus-nucleus collisions. An essential 
and here particularly relevant aspect of deconfinement is that the constituent partons of 
a deconfined medium are no longer constrained to the distribution functions of individual 
hadrons, as determined in deep inelastic scattering. The hadronic gluon distribution 
functions are strongly suppressed at high gluon momenta, with xg{x) ~ {l—xY for x — >■ 1, 
where a ~ 3 — 5 and x = kg/k^ is the relative fraction of the hadron momentum carried 
by the gluon. Removing this constraint will generally lead to harder gluons. Since J/ip 
dissociation requires hard gluons, the inelastic J/ip-hadron cross section becomes very 
small at low collision energies. Hence significant J/tp suppression in nuclear collisions 
requires the produced environment to be deconfined. 

In principle, the predicted threshold suppression of the J/ ^-hadron dissociation cross 
section can be measured directly. However, until such experiments are carried out, and 
to check the universality of the phenomenon, it is of interest to consider other processes 
in which the partonic aspects of the J/ip-hadion interactions play a decisive role. The 
photoproduction of heavy quarkonium states is, as we shall see, an excellent case at hand, 
with a quite extensive set of experimental data available over a wide region of incident 
energies. 

Within the conventional VMD approach, we can relate the reactions —>■ ipN and 
ipN — i> ipN and express the differential cross section of the forward J / -i/^-photoproduction 
on nucleons as [ITI 



{s,t = 0) = [u] — (g.^ = 0) (1) 

dt am^ \k.yN / dt 

where /c^^ = [s — {nia + mb)^][s — {rua — m6)^]/4s denotes the squared center of mass 
momentum of the corresponding reaction and F stands for the partial J/^/'-decay width; 
to simplify formulae, we shall in most expressions abbreviate J/ip by ip. The differential 
J/ip — N cross section is given by 

da^N^^N ^ I 1 iM,r.{s,t = 0)p (2) 

dt 647rm^(A^ — m^) 

with Aijpj^f denoting the invariant J/ip — N scattering amplitude; further, A = (pK/m^) 
is the nucleon energy in the quarkonium rest frame and p, K, q are the four-momenta of 
target nucleon, J/ip and initial photon, respectively. From the optical theorem in the 
form 

<N = , (3) 



Jot _ M^N 



we then obtain the well-known relation 

da^N^i^N , . _ rv. _ 3F(?/; ^ e+e~) [s - {tun + m^Y] [s - {tun - m^Y]{l + p^) / n2 
dt ^'^^-^>- 167ramv,(s-m^2)2 V'^n) ■ 

(4) 

where p is the ratio of real and imaginary parts of forward J/ip — N scattering amplitude. 
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Figure 1: Forward J/ip photoproduction data compared to our results with (solid line) 
and without (dashed line) the real part of the amplitude. The curves were obtained using 
a scahng PDF g 



By use of the operator product expansion 0-0], the J/ip — N scattering amplitude 
in the unphysical region A ~ is obtained in terms of the nucleonic gluon distribution 
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(5) 



The quantities tq = 4/(3as'"^c) and eo = mc{3as/4:Y correspond to the 'Bohr' radius and 
the 'Rydberg' energy of the lowest cc bound state J/ip, with rric for the mass of the charm 
quark. The gluon distribution g{x, e^) is renormalized at the quarkonium binding energy 
scale eo. Eq. (H) includes target mass corrections neglecting these would lead back 
to the form used in |Q. This form encounters problems in calculating the real part and 
hence would here give an incorrect threshold behavior. 

Taking into account the obvious analytical properties of the amplitude (^, one now 
relates the physical to the unphysical regions in A; together with the optical theorem (R|) 



3 



> 
o 

e 



o 

II 



T3 



7 N — > J/y N 



10 



10 r 




10 



□ CAMERINI 75 (SLAC) 
o GITTELMAN 75 (CORNELL) 
e ALBERT 79 
V BINKLEY82 
«. AID 96 (HERA HI) 
4 BREITWEG 97 (HERA ZEUS) 



10 



W, GeV 

Figure 2: The same as in Fig.l, but with the curves obtained using the PDF MRS H. 



this leads to sum rules for the J/ip — N cross section, 
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where y = niN/X and I{n) = {tt^/^ /2){32/3f[T{n + 5/2)/T{n + 5)]r3eo(mjv/eo)""^ In a 
first iteration, the solution of these sum rules can be written as a convolution of the gluon 
distribution function and the gluon- J/^/; cross section (see [0]), 
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(7) 



neglecting terms of order rrifyf'^ /eo"^. Including higher order terms, the full solution can be 
obtained iteratively, each step providing a contribution from the corresponding term of 
the hypergeometric series. Effectively, these target mass corrections change the x- variable 
in the convolution (|^) and thus the resulting threshold behaviour. 



For a given gluon distribution g(x,eo), we thus obtain and hence the imaginary 
part of the J/tp — N scattering amplitude. To determine the forward J/'?/'-photoproduction 
cross section (see Eq's. (|I]) and (0)), we need the real part of the amplitude as well. The 
high energy behaviour of the corresponding amplitude makes it possible to express this 
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Figure 3: Data for the elastic J/V'-photoproduction cross section, compared to the pre- 
diction obtained with the PDF MRS H. 

in terms of dispersion integrals with one subtraction, performed e.g. at A = @, |^: 

2A2 /-oo dX' ImM,,N{y) 



A'2 - A2 



The subtraction constant Ai^ Af(0), which is needed to estimate the behaviour of the real 
part of the amplitude near threshold is obtained by calculating the corresponding limit 
(A M^) in Eq. (j^). This is a self-consistent approach if we are given the analytical 
expression Eq. (j^) for the amplitude. A more phenomenological and perhaps also more 
realistic procedure would be to use the low-energy theorems to estimate the subtraction 
constant, following some recent work 

To complete our calculation, we have to specify the gluon distribution function of the 
nucleon and fix the overall normalization in terms of the different constants in Eq. (^) 
and related quantities. For g{x,ml), we use two parametrizations specified in |^ and 0. 
The MRS H of parametrization, in particular, takes into account the recent HERA 
results at small x and thus seems preferable for our analysis. Although all quantities in 
our formulae are in a sense 'physical' constants, uncertainties enter through the c-quark 
mass rric and the resulting values of as and eo; an addition, both the VMD model and the 
Coulombic description of the J/ ip may require corrections. We therefore treat the overall 
normalization of our result as an open constant, as was also done in 0. 

We now want to compare our results with the data available from low energy Cornell- 
SLAC up to recent high energy HERA studies These experiments measure 
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the photoproduction cross section as function of the invariant momentum transfer t in 
the 'diffractive' region at small t. One can now fit the t-dependence in the customary 
exponential form exp(— 6t) and then either extrapolate it to the unmeasurable limit 
t = at which our prediction (|1]) holds, or integrate over t from to = tmax(s) to infinity 
to obtain the 'elastic' photoproduction cross section cr^Tv^j/i/'Af (elastic). 

Figs. 1 gives the c.m.s. energy dependence of the forward differential cross section 
for the scaling PDF-parameterization g{x) = 2.5(1 — x)^ used in we show both our 
complete result and the form obtained by neglecting the real part, as in [Q. It is seen that 
the inclusion of the real part greatly improves the threshold behaviour; the agreement 
is now quite good, except for the high energy data, for which the small x behaviour of 
the PDF becomes important. Since the PDF form used here does not include this, it is 
clear that there will be deviations at high W. In Fig. 2, we then show the corresponding 
result for the new PDF MRS H, which does include the small x results from Hera. While 
the qualitative agreement is reasonable, there are definite deviations; these would become 
weaker for a less singular small x form of the PDF. Moreover, there is some experimental 
discrepancy between t-dependence of the SLAG data compared to other data in the same 
energy region; we shall return to this shortly. It is also not clear if a fine-tuning of rric 
and as would improve the situation. 




Figure 4: The energy dependence of the slope parameter for differential elastic J/ip- 
photoproduction. The solid line corresponds to a best fit with b = qq + ailnW^, giving 
Qq = —1.64 ± 0.26 and ai = 0.83 ± 0.06; for data references, see fill 



In Fig. 3, we then show our results for the exclusive J/V' photoproduction cross section 
as a function of W, compared to relevant experimental data; the theoretical error range 
is explained below. The curve was obtained for the MRS H parameterization of the PDF, 
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using the approach mentioned above. Thus, we integrate the forward differential cross 
section in the diffractive peak up to t^ax [0 with the slope parameter b as logarithmic 
function of W ('diffraction cone shrinkage'). This integrated result is much less sensitive 
to the parameter values involved and provides a good interpolation between the high 
energy (small x) data and the near-threshold behaviour of the cross section. 

Finally, we want to return shortly to the slope parameter parametrisation we have 
used above. In Fig. 4 we show the result of fitting the form 6 = oq + cti log W"^ to all 
available measurements; a best fit gives oq = —1.64 ± 0.26 and ai = 0.83 ± 0.06. It is 
seen that the SLAC data give a slope parameter considerably higher than that of the 
other experiments in the same energy region. This leads directly to the comparatively 
high values of da/dt(t = 0), which are obtained by an extrapolation to t = using 
the measured slope parameter. - The energy variation of b indicates the presence of 
contributions from 'soft' as well as from a 'hard' Pomeron Here the role of the low 



energy points is crucial; a restriction to only high energy points |T3[ could lead to an 
energy-independent slope as expected for a 'hard' Pomeron. 

In conclusion, we note that our analysis of photoproduction confirms the relation 
between the energy dependence of the cross section and the x-dependence of the gluon 
distribution function of the nucleon This relation was derived for J/'?/'-hadron in- 
teractions and enters photoproduction through VMD. The success of such a consistent 
treatment of these two reactions does not support expectations about inherently different 
behaviour |T3|] . 
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